On the dynamic universality class of the 4-state Potts model 
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We investigate the short-time critical dynamics of the 4-state Potts (FSP) model, Baxter- Wu 
(BW) model and Ising model with three-spin interactions (IMTSI) in one direction, that belong to 
the same universality class at equilibrium. The scaling behavior of the global persistence probability 
P(t) is studied in order to estimate the global persistence exponent 9 g of the models. We also 
calculate the exponent xo, that represents the anomalous dimension of the initial. Although the 
exponents 6 g and xo for the FSP model are in complete agreement with the results of the IMTSI 
model, they are completely different of that for the BW model, showing that the BW model does 
not belong to the same dynamic universality class of the FSP model. 
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I. INTRODUCTION 

Since the works by Janssen, Shaub and Schmittmann 
PJ, and Huse Q, the critical properties of statistical 
systems have been a subject of considerable interest in 
nonequilibrium physics. By using renormalization group 
methods and numerical calculations, respectively, they 
showed that there is universality and scaling behavior 
even at the early stage of the time evolution after quench- 
ing from high temperatures to the critical one. 

The dynamic scaling relation obtained by Janssen et 
al. for the kth moment of the magnetization, extended 
to systems of finite size is written as 

M<» (t, r, L, mo) = b k ^ v M {k ^ (b~ z , b l ' v r, b^L, 6 x °m ), 

(1) 

where t is the time evolution, b is an arbitrary spatial 
scaling factor, t = (T — T c )/T c is the reduced tempera- 
ture and L is the linear size of the lattice. The exponents 
(3 and v are as usually the equilibrium critical exponents 
associated respectively with the order parameter and the 
correlation length, and z is the dynamical exponent char- 
acterizing time correlations in equilibrium. 

Besides to avoid the well-known problem of the "criti- 
cal slowing down" , characteristic of the equilibrium, and 
to provide an alternative way to obtain the familiar set of 
static critical exponents and dynamic critical exponent z, 
in this new universal regime a new dynamic critical ex- 
ponent 9 is found by following the above scaling law for 
the order parameter at the critical temperature (r = 0) 

M(t) ~ m t e , (2) 

where mo represents the initial magnetization of the sys- 
tem. This new index, independent of the previously 
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known exponents characterizes the so-called "critical ini- 
tial slip" , the anomalous increase of the order parameter 
when the system is quenched to the critical temperature 
T c . In addition, a new critical exponent xq, which rep- 
resents the anomalous dimension of the initial magneti- 
zation is introduced to describe the dependence of the 
scaling behavior on the initial conditions. This exponent 
is related to 9 as 

x = Oz + (3 jv. (3) 

As proposed by Majumdar et al. another dynamic 
critical exponent can be obtained in the study of sys- 
tems far from equilibrium. By studying the behavior of 
the global persistence probability P(t) that the order pa- 
rameter has not changed its sign up to time t, they shown 
that P(t) should behave, at the critical temperature, as 

P(t)~f-*«, (4) 

where 9 g is the global persistence exponent. They also 
argued that, if the time evolution would be a Markovian 
process, then the exponent 9 g should obey the equation 

1 

9 g z = -9z + --^-. (5) 

Z V 

However as shown in several works 0, IE El ESS El El 
lialiallilla the exponent 9 g is an independent critical 
index closely related to the non-Markovian process. 

In this work, we perform short-time Monte Carlo sim- 
ulations to explore the scaling behavior of th e gl obal per- 
sistenceprobability Pit) for the FSP model E3 , BW 
model [HEl and IMWTSI [H E3, that belong to the 
same universality class at equilibrium. We also calculate 
the exponent xq of these models through the Eq. 10) • 
However, before calculating this exponent, we decide to 
reobtain more reliable estimates of the exponents 9 and 
z for the FSP and IMTSI models. 

The paper is organized as follows. In the next section 
we present the models. In Section ITTT1 we show the short- 
time scaling relations and present our results. Finally, in 
Section llVl we present our conclusions. 
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II. THE MODELS 

The q— state Potts model which is a simple extension 
of the Ising model, has a rich phase diagram with first 
order phase transitions when q > 4 and second order 
phase transitions when q < 4. Its Hamiltonian is given 
by 

pH = -K^8<, i( , v (6) 

{id) 

where (3 — l/fc^T and ks is the Boltzmann constant, 
represents nearest-neighbor pairs of lattice sites, 
K is the dimcnsionless ferromagnetic coupling constant 
and <ii is the spin variable which takes the values di — 
0, • • • , q — 1 on the lattice site i. It is well known that 
the critical point of this model is given by 01 

K c = log(l + V?), (?) 
and its order parameter is defined as 

M(f) = L^T) \ I>W - !) } ( 8 ) 

where L is the linear size of the lattice and d is the di- 
mension of the system. 

The BW model is defined by the Hamiltonian 

-PH = K J2 a i a j a k, (9) 

(i,3,k) 

where <Xj = ±1 are Ising spin variables located at the sites 
of the triangular lattice and the sum is over all elementary 
triangles. 

The Hamiltonian of the Ising model with three-spin 
interactions in the horizontal direction and simple two- 
spin interactions in the vertical direction is given by 

— f3H = ^ {K h (Tij<T t+1 .j<T i+ 2,j + K v ai^a lt j + i\ , (10) 

where the sum is over all sites of a square lattice, Kh 
and K v are the coupling constant in the horizontal and 
vertical directions, respectively, and a^j = ±1 are Ising 
spin variables located at the sites of the lattice. 

Both BW and IMTSI (for the isotropic case, K h = 
K v = K) models undergo a continuous phase transition 
at critical temperature K c — 0.51n(v2 — 1) which is the 
same temperature as the Ising model on a square lattice. 
The order parameter of these models is defined as, 



M(t) = 




The BW and IMTSI models present semi-global up- 
down spin reversal symmetry, i.e., their Hamiltonians are 
invariant under reversal of all the spins belonging to two 



of three sublattices into which the original lattice can be 
decomposed. 

The ground state of these models are fourfold degen- 
erate, being that the possible spin configurations of the 
BW and IMTSI models consist of repetitions of the pat- 
terns {+, +, +}, {+, -, -}, {-, +, -} or {-, -, +}. The 
main difference between them is that the BW model is 
defined on a triangular lattice whereas in the FSP and 
IMTSI models the spins are located on a square one. 

From the degeneracy and symmetry considerations, 
one can conjecture that both models belong to the same 
universality class at equilibrium as the FSP model with 
exponents [22| 

12 1 
/?=— , u = a=-, and J? = — . (12) 

However, when deepenly study one these models, dif- 
ferences among sub-dominant exponents appear. These 
exponents are responsible for scale corrections, which can 
be related to the character of the convergence of the re- 
sults obtained through the finite size scaling. These dif- 
ferences arised on the first time in the study of the con- 
forme invariance. 

In addition, when the dynamic critical behavior of 
these models is take into account, the dynamic expo- 
nents found for the BW model are completely differ- 
ent from that values of the FSP model universality class 

III. RESULTS AND DISCUSSIONS 

In our Monte Carlo simulations, we consider two- 
dimensional lattices with periodic boundary conditions. 
The dynamic evolution of the spins in local and updated 
by the heat-bath algorithm at the critical point K c . In 
order to check the possibility of finite-size effects, we con- 
sider three different lattice sizes, L = 120, 180 and 240 
and the exponents were obtained from five independent 
bins of 20000 samples each one. 



A. Global persistence exponent 6 g 

The global persistence probability P(t) can be defined 

as 

P(t) =1-J>(f) (13) 
t'=i 

where p{t') is the fraction of the samples that have 
changed their state for the first time at the instant t' . 
The dynamical exponent g that governs the behavior 
of P(t) at criticality is obtained through the power law 
behavior given by 

P(t)~t~ e >. (14) 
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In order to obtain the exponent 9 g , the initial con- 
figuration of the system should be carefully prepared 
with a precise and small value of too. After estimat- 
ing 9 g for a number of too values, its final value is ob- 
tained from the limit too — * 0. In this work, we used 
4 • 1(T 4 < too < 5 • 1(T 3 . 

As we are consider three models and two different or- 
der parameters, it is worth to explain how obtain too in 
each case. For the FSP model, each site on the lattice is 
occupied by a spin variable which takes the values a = 0, 
1, 2 or 3 with equal probability. After, the magnetization 
is measured by using the Eq. JSJ) and then, the variables 
in the sites are randomly chosen up to obtain a null value 
for the magnetization. The last procedure is to change 
8 sites on the lattice in order to obtain the initial mag- 
netization desired. When using the Eq. JjjJ, the initial 
magnetization is given by 



toq 



4(5 



(15) 



and a value of to is obtained choosing 8 sites occupied 
by a = 0, 2 or 3 and substituting them by a = 1. For 
the Eq. Ullll . toq is given simply by 



(16) 



and to obtain a value of Too, we choose randomly 8/2 
sites occupied by a = — 1 and substitute them by a = 1. 

In Fig. we show the behavior of the global persis- 
tence probability for L = 240 and a small value of Too 
for the FSP (on top), BW (on middle), and IMTSI (on 
bottom) models, in double-log scales. The error bars, 
calculated over five sets of 20 000 samples are smaller 
than the symbols. 

The inset in Fig. ^ display the estimates of 9 g for 
different values of mo and the limit procedure too — » 
for the models. 

In Table [I] we show the extrapolated values of 6 g for 
180 and 240 for the FSP, BW and IMTSI mod- 



L = 120, 
els. 



TABLE I: The global persistence exponent 6 g from the power 
law behavior for the FSP, BW, and IMTSI models. 



Models 


L = 120 


L = 180 


L = 240 


FSP model 


0.469(4) 


0.472(6) 


0.475(5) 


BW model 


0.620(5) 


0.618(5) 


0.619(4) 


IMTSI model 


0.469(6) 


0.470(5) 


0.471(5) 





0.50 r 
0.4S ; 

0.46 - 
0.44 - 
0.42 - 
0.40 - 



= 0.440(8) 



0.001 0.002 0.003 0.004 0.005 



= 0.619(4) 



0.001 0.002 0.003 0.004 0.005 



0.001 0.002 0.003 0.004 0.005 

m o 



L = 240 



m B = 0.000625 



BW model 



L = 240 

m =0.0003125 




L = 240 

m =0.0003125 



The values found for the three different lattice sizes 
are shown in Tableland indicate that the finite size ef- 
fects are less than the statistical errors. The discrepancy 
among the results of the BW model and the results of the 
FSP and IMTSI models is noteworthy. Discrepancies be- 
tween the BW and FSP models, were also observed by 
Arashiro and Drugowich de Felicio [23j for the dynamical 
exponent 8, by Hadjiagapiou et al. 24] for the exponents 



FIG. 1: The time evolution of the global persistence proba- 
bility P(t) for L = 240 for the FSP (on top), Bw (on middle), 
and IMTSI(on bottom) models. The error bars calculated 
over 5 sets of 20 000 samples are smaller than the symbols. 
The inset in each figure shows the exponent 6 g for different 
initial magnetizations, as well as its extrapolated value. 



4 



9 and z, and by Chatelain [25( for the exponent X/z and 
for the asymptotic value of the fluctuation-dissipation ra- 
tio Xrr,- 



B. Dynamic critical exponents 9 and z 

As conjectured by Janssen et al. on basis of renor- 
malization group techniques and by Huse @] through 
numerical calculations, in the short-time regime the or- 
der parameter decays algebraically as shown in the Eq. 
(J2J, where 9 is the dynamic critical exponent that char- 
acterizes the anomalous behavior of the magnetization. 
Formely, a positive value was always associated to this 
exponent HE IHIH |H H3, IH IH and the phenomenon 
was known as critical initial increase. However, as shown 
in some papers, there are models in which the exponent 
9 can have a negative value, for instance, the tricritical 
Ising model, analitically studied by Janssen and Oerding 
|33|| and numerically confirmed by da Silva et al. |34| . 

As shown in several works, this is the case of the FSP 
HH, BW IHH, and IMTSI jH models. 

Although the estimates of the exponents 9 and z for 
the BW model seem to be with a good precision, the re- 
sults for the IMTSI model seem to have a large error. In 
addition, the estimates obtained for 9 in Ref. [3^ show 
a considerable difference between the two techniques em- 
ployed to the FSP model. So, in order to estimate a more 
precise value of the exponent xq, we reobtain the expo- 
nents 9 and z for the FSP and IMTSI models by using 
the time correlation of the magnetization p9| 



C(t) = (M(O)M(i)) ~ t e 
and mixed initial conditions |37| 

(AP(t)) 



(M(t)) 



/m =Q 
2 

mo=l 



t d/z 



(17) 



(18) 



respectively. In Eq. (|17fl . the average is taken over a 
set of random initial configurations. Initially, this ap- 
proach had shown to be valid only for models which ex- 
hibit up-down symmetry [30j. Nevertheless, it has been 
demonstrated recently that this approach is mo re g eneral 
and can include models with other symmetries 38] . This 
method has several advantages when compared to other 
approachs, for instance, the exponent 9 can be directly 
calculated, without the need of careful preparation of the 
initial states nor of the limiting procedure [see Eq. @], 
the only requirement being that (M(0)) = 0. 

In Fig. |21 we show the time dependence of the time 
correlation C(t) in double-log scale for the FSP (on top) 
and IMTSI (on bottom) models, respectively, for L = 
240. 

The slope of these curves is shown in Table ITT1 as well 
as the estimates for L = 120 and L = 180. 

On the other hand, the dynamical exponent z was ob- 
tained by combining results from samples submitted to 
different initial conditions (see Eq. Q18[). where d = 2 is 



FSP model 



L = 240 

6 = - 0.046(9) 




L = 240 

9 = - 0.047(8) 



FIG. 2: The time correlation of the order parameter on log 
scales, for the FSP (on top) and IMTSI (on bottom) models. 
Error bars were calculated over five sets of 20 000 samples. 



TABLE II: The exponent for the FSP and IMTSI models. 



L 


FSP model 


IMTSI model 


120 


-0.046(8) 


-0.047(7) 


180 


-0.047(8) 


-0.046(7) 


240 


-0.046(9) 


-0.047(8) 



de dimension of the system). This approach has proved 
to be very efficient in estimating the exponent z for sev- 
eral models |23|, |32j, 1341 1371 |39j ■ The time evolution of 
is shown on log scales in Fig. [3] for L = 240 for the FSP 
(on top) and IMTSI (on bottom) models. 

Taking into account the values of the ratio d/z, es- 
timated from the slope of these curves, the exponent z 
can be easily found. Our estimates of this exponent for 
the FSP and IMTSI models are shown in Table |TTT| for 
L = 120, 180, and 240. 

As shown in Tables ITTl and IIIII the estimates of the 



FSP model 




0.871(2) 




FIG. 3: The time evolution of F 2 (t) for L = 240 for the 
FSP (on top) and IMTSI (on bottom) models. Each point 
represents an average over five sets of 20 000 samples and the 
error bars are obtained of them. 



TABLE III: The exponent z for the FSP and IMTSI models. 



L 


FSP model 


IMTSI model 


120 


2.294(7) 


2.293(5) 


180 


2.294(5) 


2.290(8) 


240 


2.296(5) 


2.292(4) 



and the estimates obtained by Hadjiagapiou el al. |24j . 

9 = -0.185(2) and z = 1.994(24). (20) 

C. Anomalous dimension xq 

Finally, we calculate the value of the anomalous di- 
mension xo of the order parameter for the FSP, BW, and 
IMTSI models. This exponent, which is introduced to 
describe the dependence of the scaling behavior of the 
initial conditions, is related to the exponents 9, z, and 
/3/Vby 



Xq = 0Z + j3 /v. 



(21) 



In Table llVl we show the estimates of xq, by using the 
values of the Tables HU and EH for L = 240 (FSP and 
IMTSI models), the results of the Eqs. (O and (j2Uj) . 
and the conjectured values of fj and v given by Eq. i|12[l . 



TABLE IV: The exponent x for the FSP, BW, IMTSI mod- 
els. 



Models 


x 


FSP model 


0.019(21) 


BW model [23J 


-0.302(6) 


BW model [24J 


-0.244(8) 


IMTSI model 


0.017(18) 



Because the value of the exponent 9, the anomalous 
dimension of the BW model is also different of the esti- 
mates for the FSP and IMTSI models. While xo is neg- 
ative for the BW model, it have a null value for the FSP 
and IMTSI models whose meaning is the presence of the 
marginal operator, i.e., the operator which have the scal- 
ing dimension equal to dimensionality of the system and 
whose effect is not modified under renormalization-group 
operations. An alike value was found recently by da Silva 
and Drugowich de Felfcio [35| for the FSP model. By us- 
ing mixed initial conditions for the magnetization, they 
found Xo = 0.0176(75) that is in very good agreement 
with our results. 



dynamical exponents 9 and z for the FSP and IMTSI 
models are in complete agreement by considering the er- 
ror bars. Already the study of the dynamical behavior 
of the BW model has shown that the exponent 9 is very 
different from that exponent of the FSP and IMTSI mod- 
els and that there are still controversies about the result 
of the exponent z, with values that extends from ~ 2 to 
~ 2.3. For instance, Arashiro and Drugowich de Felfcio 
have obtained p3|. 

9 = -0.186(2) and z = 2.294(6) (19) 



IV. CONCLUSION 

We estimate the dynamic exponent 9 g of the FSP, BW 
and IMTSI models using the time-evolution of the global 
persistence probability that the magnetization has not 
changed its sign up to time t. The value of 9 g found for 
the BW model is completely different of that found for 
the FSP and IMTSI models. On the other hand, our re- 
sults obtained for the FSP and IMTSI models are in good 
agreement each other. We also reobtained the dynamic 
exponents 9 and z of the FSP and IMTSI models in order 
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to calculate de exponent xq. Our estimates of these ex- 
ponents confirm that the FSP and IMTSI models belong 
to the same universality class far from equilibrium. Fol- 
lowing, we also calculated the exponent x of the FSP, 
BW, IMTSI models. Our estimates of x for the FSP 
and IMTSI models are very closed to zero, whereas for 
the BW model it is very negative. These results, together 
the estimates of the exponents g , 6, and z allow us assert 
that, although the BW model belong to the universality 
class of the 4-state Potts model at equilibrium, it seems 
not be the case far from equilibrium. These results sug- 



gest that the short time simulations get to reveal some 
properties that in the equilibrium stage are very difficult 
of note. 
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